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Abstract 

We discuss the properties of an ideal relativistic gas of events possessing 
Bose-Einstein statistics. We find that the mass spectrum of such a system is 
bounded by fi < m < 2M//j,k, where fi is the usual chemical potential, M is 
an intrinsic dimensional scale parameter for the motion of an event in space- 
time, and hk is an additional mass potential of the ensemble. For the system 
including both particles and antiparticles, with nonzero chemical potential fx, 
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the mass spectrum is shown to be bounded by \fi\ < m < 2M//j,k, and a special 
type of high-temperature Bose-Einstein condensation can occur. We study this 
Bose-Einstein condensation, and show that it corresponds to a phase transition 
from the sector of continuous relativistic mass distributions to a sector in which 
the boson mass distribution becomes sharp at a definite mass M/^ik- This 
phenomenon provides a mechanism for the mass distribution of the particles to 
be sharp at some definite value. 

Key words: special relativity, relativistic Bose-Einstein condensation, mass distribu- 
tion, mass shell 
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1 Introduction 

There have been a number of papers in the past H, ||, ||, which discuss the properties 
of an ideal relativistic Bose gas with nonzero chemical potential fi. Particular attention 
has been given to the behavior of the Bose-Einstein condensation and the nature of 
the phase transition in d space dimensions 0, ^ . The basic work was done many years 
ago by Jiittner ||, Glaser J7|, and more recently by Landsberg and Dunning-Davies 
H and Nieto ||. These works were all done in the framework of the usual on-shell 
relativistic statistical mechanics. 

To describe an ideal Bose gas in the grand canonical ensemble, the usual expression 
for the number of bosons N in relativistic statistical mechanics is 



where V is the system's three-volume, E k = y/k 2 + m 2 and T is the absolute tem- 
perature (we use the system of units in which h = c = ks = I; we also use the 
metric g^ u = (—,+,+, +)), and one must require that fi < m in order to ensure a 
positive-definite value for n k , the number of bosons with momentum k. Here N is 
assumed to be a conserved quantity, so that it makes sense to talk of a box of iV 
bosons. This can no longer be true once T ~ m [10]; at such temperatures quantum 
field theory requires consideration of particle-antiparticle pair production. If N is the 
number of antiparticles, then N and N by themselves are not conserved but N — N 
is. Therefore, the high-temperature limit of (1.1) is not relevant in realistic physical 
systems. 

The introduction of antiparticles into the theory in a systematic way was made 
by Haber and Weldon [1(| |j~If . They considered an ideal Bose gas with a con- 
served quantum number (referred to as "charge") Q, which corresponds to a quan- 
tum mechanical particle number operator commuting with the Hamiltonian H. All 



N = VY j n k = VY J 
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thermodynamic quantities may be then obtained from the grand partition function 



Tr {exp [—{H — pQ)/T]} considered as a function of T, V, and p, |12| . The formula 
for the conserved net charge, which replaces (1.1), readsQ [10] 
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In such a formulation a boson-antiboson system is described by only one chemical 
potential /i; the sign of p indicates whether particles outnumber antiparticles or vice 
versa. The requirement that both and fik be positive definite leads to the important 
relation 

| /j- | < m. (1.3) 

The sum over k in (1.2) can be replaced by an integral, so that the charge density 
p = Q/V becomes 
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which is an implicit formula for p as a function of p and T, and in the region T » m 
reduces to 

(1-5) 

For T above some critical temperature T c , one can always find a pi < m) such 
that (1.4) holds. Below T c no such /i can be found, and (1.4) should be interpreted 
as the charge density of the excited states: p — p , where p is the charge density 
of the ground state [10] (with k = 0; clearly, this state is given with zero weight in 
the integral (1.4)). The critical temperature T c at which Bose-Einstein condensation 
occurs corresponds to pi = ±m (depending on the sign of p.) Thus, one sets \p,\ = m 
in (1.4) and obtains, via (1.5) (provided that |p| >> m 3 ), 



Tc= m. (1 . 6) 

V m 

Below T c , (1.4) is an equation for p — p Q , so that the charge density in the ground 
state is 

Po =p[l-(T/T c ) 2 ]. (1.7) 

It follows from Eq. (1.6) that any ideal Bose gas will condense at a relativistic 
temperature (T c >> m), provided that |p| >> m 3 . 

Recently the analogous phenomenon has been studied in relativistic quantum field 
theory [[11], [13|, [TJ], [H| . For relativistic fields Bose-Einstein condensation occurs at high 



-'^One uses the standard recipe according to which all additive thermodynamic quantities are 
reversed for antiparticles. 



3 



temperatures and can be interpreted in terms of a spontaneous symmetry breaking 
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The extension and generalization of Bose-Einstein condensation to curved space- 
times and space-times with boundaries has also been the subject of much study. The 
non-relativistic Bose gas in the Einstein static universe was treated in ref. |]]. The 
generalization to relativistic scalar fields was given in refs. [jl6| , [17|j . The extension to 
higher dimensional spheres was given in ref. [18|. Bose-Einstein condensation on hy- 



perbolic manifolds |19| |, and in the Taub universe 0] has also been considered. More 
recently, by calculating the high-temperature expansion of the thermodynamic poten- 
tial when the boundaries are present, Kirsten [21] examined Bose-Einstein condensa- 
tion in certain cases. Later work of Toms [2^] showed how to interpret Bose-Einstein 
condensation in terms of symmetry breaking, in the manner of flat space-time cal- 
culations [1TTL |13f. The most recent study by Lee et al. [p3l showed how interacting 



scalar fields can be treated. Bose-Einstein condensation for self-interacting complex 
scalar fields was considered in ref. 
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In the present paper we consider a relativistic Bose gas within the framework of a 
manifestly covariant relativistic statistical mechanics ]55, 5fj, |27j. We shall review this 
framework briefly in the next section. We obtain the expressions for characteristic 
thermodynamic quantities and show that they coincide with those of the relativistic 
on-shell theory, except for the value of the average energy (which differs by a factor 
2/3). We introduce antiparticles and discuss the properties of Bose-Einstein conden- 
sation in such a particle-antiparticle system. We show that it corresponds to a phase 
transition to a high-temperature form of the usual on-shell relativistic kinetic theory. 



2 Relativistic TV-body system 

In the framework of a manifestly covariant relativistic statistical mechanics, the dy- 
namical evolution of a system of N particles, for the classical case, is governed by 
equations of motion that are of the form of Hamilton equations for the motion of N 
events which generate the space-time trajectories (particle world lines) as functions 
of a continuous Poincare-invariant parameter r, called the "historical time" [2S[ |2S| . 
These events are characterized by their positions = (t, q) and energy-momenta 
= (E, p) in an 8iV-dimensional phase-space. For the quantum case, the system 
is characterized by the wave function Vv(?i> <?2, • • • , <1n) G L 2 {R m ), with the measure 
(ftqidfiqz ■ ■ ■ d A q N = d m q, (qt = gf ; fi = 0, 1, 2, 3; i — 1, 2, . . . , N), describing the dis- 
tribution of events, which evolves with a generalized Schrodinger equation [29|| . The 
collection of events (called "concatenation" J3(|) along each world line corresponds 
to a particle, and hence, the evolution of the state of the iV-event system describes, 
a posteriori, the history in space and time of an iV-particle system. 
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For a system of N interacting events (and hence, particles) one takes |29 



* = £^ + V(fc,ft,...,?tf), (2.1) 

i 

where M is a given fixed parameter (an intrinsic property of the particles), with the 
dimension of mass, taken to be the same for all the particles of the system. The 
Hamilton equations are 

dq- = dK_ = 

(It dp ifi M' 

dp? dK dV 
dr dq^ dq^ ' 
In the quantum theory, the generalized Schrodinger equation 



(2.2) 



d 

i-Q-.AiQi, 92, • • • , ?tv) = Kij) T (qx, q 2 , ■ ■ ■ , q N ) (2.3) 

describes the evolution of the iV-body wave function ip T (qi, q 2 , ■ ■ ■ , ?v)- To illustrate 
the meaning of this wave function, consider the case of a single free event. In this 
case (2.3) has the formal solution 

Mq) = (e- lKoT ^o)(q) (2.4) 

for the evolution of the free wave packet. Let us represent i/j T (q) by its Fourier 
transform, in the energy-momentum space: 



MQ) = J^y 2 J d*pe-^e^Mp), (2-5) 

where p 2 = p^p^, p-q = p^q^, and ipo(p) corresponds to the initial state. Applying the 
Ehrenfest arguments of stationary phase to obtain the principal contribution to ip T {q) 
for a wave packet at one finds (p% is the peak value in the distribution ipoip)) 

q» * §r, (2.6) 

consistent with the classical equations (2.2). Therefore, the central peak of the wave 
packet moves along the classical trajectory of an event, i.e., the classical world line. 



In the case that p° c = E c < 0, we see, as in Stueckelberg's classical example f28 
that 

dt c E 

— — — -rj < 0. 

dr M 



It has been shown ]30| in the analysis of an evolution operator with minimal electro- 
magnetic interaction, of the form 

(p-eA(q)) 2 
2M 
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that the CPT-conjugate wave function is given by 

W CPT (t,q)=VvH,-q), (2.7) 
with e — > —e. For the free wave packet, one has 

tf PT (q) = j^y 2 f d'pe-^e-^Mp)- (2-8) 
The Ehrenfest motion in this case is 

if E c < 0, we see that the motion of the event in the CPT-conjugate state is in the 
positive direction of time, i.e., 



dt r E r E r 



dr " M M 



(2.9) 



and one obtains the representation of a positive energy generic event with the opposite 
sign of charge, i.e., the antiparticle. 



2.1 Ideal relativistic Bose gas 

To describe an ideal gas of events obeying Bose-Einstein statistics in the grand canon- 



ical ensemble, we use the expression for the number of events found in |25 
N = V® V n k » = V - 2 , 



(2.10) 



where is the system's four- volume, fix the additional mass potential [25|, which 
we shall take, in order to simplify subsequent consideration, to be a fixed parameter 
(which determines an upper bound of the mass distribution in the ensemble we are 
studying, as we shall see below), and m? = —k 2 = —k^k^. To ensure a positive- 
definite value for n^, the number of bosons with four-momentum k^, we require 
that 

(2.11) 



m 



m — n — hk 



2M 



> 0. 



The discriminant for the l.h.s. of the inequality must be nonnegative, i.e., 

M 



2fi K 



(2.12) 



For such /i, (2.11) has the solution 




< m < 



M 



2/iyU 
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(2.13) 
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For small hhk/M, the region (2.13) may be approximated by 

fx<m< — . (2.14) 

One sees that plays a fundamental role in determining an upper bound of the 
mass spectrum, in addition to the usual lower bound m > //. In fact, small fix admits 
a very large range of off-shell mass, and hence can be associated with the presence of 



strong interactions |31 



Replacing the sum over (2.10) by an integral, one obtains for the density of 
events per unit space-time volume n = N/V^ p2| , 



„ = JL r dm r — '" a sinh f — , (2.i5) 

47T 3 J mi J-oo e (mcoBb.p-vi-it K %z)/T _ I 

where mi and m 2 are defined in Eq. (2.13), and we have used the parametrization 

m 

p° = m cosh /3, 

p 1 = m sinh (3 sin 9 cos 4>, 

p 2 = m sinh (3 sin 9 sin <p, 

p 3 = m sinh (3 cos 0, 

< 6> < 7T, < < 27T, -oo < (3 < oo. 

In this paper we shall restrict ourselves to the case of high temperature alone: 

M 

T»—. (2.16) 

fJ>K 

It is then possible (for simplicity here; see footnote 3 on Eq. (3.9)) to neglect indis- 
tinguishability of bosons in the integrand^ and to rewrite (2.15) in the form 



e /VT 
47T 3 



/ m 3 dm smh 2 p dPe- mcoah P /T e flKm /2MT , (2.17) 



which reduces, upon integrating out /?, to p7 



n 



471-3 



P dm m 2 ^ W eWtmVaMT (2 .18) 

imi \1 J 



where K v (z) is the Bessel function of the third kind (imaginary argument). Since 

m < m 2 < 2M/ 

VKm 2 < fi K m/fi K ) 2 = ™_ 

2MT ~ 2MT Tfi K ' 1 ' J 



2 The limits of integration remain as a vestige of the Bose-Einstein distribution; they play a central 
role in the results, as we shall see below. 
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u m 2M , 

t^t^Wk <<x - (2 ' 20) 



in view of (2.16), and also 
Therefore, one can neglect the exponentials in Eq. (2.18), and for Ki(m/T) use the 

K v (z)~±.r(v)(~) \ z«l. (2.21) 



asymptotic formula |33 



2 



Thus, we obtain 



4tt 3 J mi 2tt 3 \fi K / f M 

From this equation, one can identify the high-temperature mass distribution for the 
system we are studying, f{m) ~ m, so that 



»4r*»-Sffl'i'-¥. (2-22) 



, e _ CI dm m e+1 2 mf 2 - mj+ 2 



In particular, 



^-©V^)- (2 - 25) 

Extracting the joint distribution for j3 and m from (2.17) in the same way, let us also 
obtain the average values of the energy and the energy squared for high T. First, 

, x , . C 2 m A dm sinh 2 (3 cosh (3d(3e- mcosh ^ T 

(E) 5 (mcoshfl S J - Cm3dmsi ; h . ^ ^_ . (2 . 26) 

Integrating out /3, one finds 

zpv ^ 1 C dm m 4 [ir 3 (m/T) - ^(m/T)] 

^ ; 4T J^dmm^K^m/T) ' 1 j 

It is seen, with the help of (2.21), that it is possible to neglect K\ in comparison with 
K% in the numerator of (2.27) and obtain, via (2.21), 

, v 1 C 2 dmm 4 K 3 (m/T) 

(E) = — p- 3 ' ~ 2T, 2.28 

4T J mi 2 dm m z Ki{m/T) 
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in agreement with refs. |23, ^6], Similarly, one obtains 

J"™ 2 m 5 dm sinh 2 /? cosh 2 f3d/3e- mcosh ^ T 



(E ) = (m 2 cosh 2 0) 



f™ 2 m 3 dm sinh 2 /3d/?e- mcosh 0/ T 



C dm[m%(m/r) + 3Tm 3 ^ 2 (m/r)] Q dmm*K 2 (m/T) ^ 2 

CI dm m 2 K x (m/T) J™ 2 dm m 2 ^ (m/T) 

(2.29) 

Let us assume that the average (p^p v ) has the form 

(p/y) = auV + bg» u , (2.30) 

where = (1, 0) in the local rest frame. The values of a and 6 can be then calculated 
as follows: for \i = v = one has (Qo ) 2 ) = a — b, while contraction of (2.30) with g^ u 
gives —g^ipuPv) — a — 46. The use of the expressions (2.29) for ((p ) 2 ) = (-E 2 ), and 
(2.25) for —g^ipfjPv) = ( m2 ) yields 

f a-b = QT 2 , 

{ a -46 = 2(^) 2 (l-^/M), 

so that 

2 /M\ 2 



a = 8T 



2 



^ , (2-31) 



2T 2 --( —) (l-Og). (2.32) 



3 \fi K J V M / 

2 

3 V/i^J V* M J 
For T >> M/hk-i it is possible to take a = 8T 2 , 6 = 2T 2 , and obtain, therefore, 

(pV) = 8T 2 uV + 2T 2 g flu . (2.33) 

To find the expressions for the pressure and energy density in our ensemble, we 
study the particle energy-momentum tensor defined by the relation |26| 



T^(q) = E / ^W^ 5 ^ ~ 9<(r))> (2 - 34) 



in which M/^k is the value around which the mass of the bosons making up the 
ensemble is distributed, i.e., it corresponds to the limiting mass-shell value when the 
inequality (2.12) becomes equality. Upon integrating over a small space-time volume 
AV and taking the ensemble average, (2.34) reduces to p6| 



Fn = -^n(p^). (2.35) 

M/fiK 
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In this formula Tav is the average passage interval in r for the events which pass 
through the small (typical) four-volume AV^ in the neighborhood of the _R 4 -point. The 
four- volume AV is the smallest that can be considered a macrovolume in representing 
the ensemble. Using the standard expression 

(T"") = (p + p)u^u u + pg» v , (2.36) 

where p and p are the particle pressure and energy density, respectively, we obtain 



Tav M 2pp K 4 
P = —3 V 1 T7~ T » P = 3p- 2.37 

To interpret these results we calculate the particle number density per unit three- 
volume. The particle four-current is given by the formula J26J 



nq) = E / dT jfh 6 * {q " %(r)) ' (2 ' 38) 



which upon integrating over a small space-time volume and taking the average reduces 
to 

(J*) = wr^ n ^y> ( 2 - 39 ) 

M/fJLK 



then 



so that 



No = (J°) = ^T^n(E), (2.40) 



Wo =^^l-^, (2.41) 
and we recover the ideal gas law 

V = N T. (2.42) 

Since, in view of (2.13), 



2M 2/^ 

A/1 tt~ = Am 

is a width of the mass distribution around the value M/p K , Eqs. (2. 37), (2. 41) can 
be rewritten as 



T AV Am 4 
P = — ^ — T , p = 3p, 
In 4 



JVo = (2.43) 

In ref. |34| we obtained a formula which relates the average passage interval in r, 
Tav, with a width of the mass distribution Am (in order that results of the off-shell 
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theory coincide with those of the standard on-shell theory in the sharp-mass-shell 
limit): 

T AV Am = 2tt. (2.44) 

One can understand this relation, up to a numerical factor, in terms of the uncertainty 
principle (rigorous in the L 2 (R 4 ) quantum theory) AE ■ At ~ 1/2. Since the time 
interval for the particle to pass the volume AV^ (this smallest macroscopic volume 
is bounded from below by the size of the wave packets) At = E/M At, and the 
dispersion of E due to the mass distribution is AE ~ mAm/E, one obtains a lower 
bound for TavAiti of order unity. 

Thus, with (2.44) holding, the formulas (2.43) finally reduce to 

P = ^2, P = 3p, (2-45) 
N = ^, (2.46) 

which are the standard expressions for high temperatures JEJ . Thus, the formulas for 
characteristic thermodynamic quantities and the equation of state for a relativistic 
gas of off-shell events coincide with those of the relativistic gas of on-shell particles, 
except for the expression for the average energy which takes the value IT in the 
relativistic gas of events, in contrast to 3T, as for the high-temperature limit of the 
usual theory |36|. Experimental measurement of average energy at high temperature 



can, therefore, affirm (or negate) the validity of the off-shell theory. There seems 
to be no empirical evidence which distinguishes between these results at the present 
time. The quantity a = M c 2 /ksT, a parameter which distinguishes the relativistic 
from the nonrelativistic regime (see, e.g., |2D|]) is very large for M of the order of 



the pion mass, at ordinary temperatures; the ultrarelativistic limit corresponding to 
a small becomes a reasonable approximation for T ~ 10 13 K. 



3 Introduction of antiparticles 



The introduction of (positive energy) antiparticles into the theory as negative energy 
events in the CPT-conjugate state leads, by application of the arguments of Haber 
and Weldon jnj, or Actor [37|, to a change in sign of \x in the distribution function 
for antiparticles. We therefore write down the following relation which represents the 



11 



analog of the formula (1.2): (] 



N 



1 



1 



(3.1) 



We require that the both n^'s in Eq. (3.1) be positive definite. In this way we obtain 
the two quadratic inequalities, 



m 



m + n — Hk 



m 
2M 

2 

m 
2M 



> 0, 

> 0, 



(3.2) 



which give the following relation representing the nonnegativeness of the correspond- 
ing discriminants: 

MM , nn . 

< a < . 3.3 

2fi K ~ ^ ~ 2fi K y 1 

It then follows that we must consider the intersection of the ranges of validity of the 
two inequalities (3.2). Indeed, if each inequality is treated separately, there would be 
some values of m for which one and not another would be physically acceptable. One 
finds the bounds of this intersection region by solving these inequalities, and obtains^ 





< m < 



which for small |/i|/i^/M reduces, as in (2.14) in the no-antiparticle case, to 

2M 



(3.4) 



< m < 



(3.5) 



Replacing summation in (3.1) by integration, we now obtain a formula for the 
number density: 



n 



471-3 



pni2 roc 

' rr? dm \ sinh 2 (3d(3 

mi J —oo 



1 



1 



e (mcosh/3- j u-/i ff 2t_)/T _ ^ e ( mcos h/3+ At - MK |i ? )/T _ ^ 

(3.6) 



3 As for the nonrelativistic theory, the "free" distribution functions describe quasiparticles in 
a form which takes interactions into account entering through the chemical potential. Since, by 
definition of a good quasiparticle, it is not frequently emitted or absorbed. We therefore consider 
the particles and antiparticles as two species. Since the particle number is determined by the 
derivative of the free energy with respect to the chemical potential, u must change sign for the 
antiparticles [0. Similarly the average mass (squared) is obtained by the derivative with respect 
to /iif ^5|; since the mass (squared) of the antiparticle is positive, nx does not change sign. 

4 This is actually the solution of one of the inequalities (the most restrictive), depending on the 
sign of ji. 
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where mi and ni2 are denned in Eq. (3.4), which for large T reduces, as above, to 



n = = T / dm m 2 ^ - e^ m ' 2MT . 

4^ J mi \T J 

Now, using the estimates (2. 19), (2. 20), and sinh(/i/T) = \ijT for jjl/T « 1, we 
obtain (in place of (2.22)) 



In (3.7) n is a conserved net event charge, the sign of fi indicating whether particles 
outnumber antiparticles or vice versa. Similarly, one obtains 

P = 2p(/i), 
P = 2 P(X>> 

where p(fi) and p(/i) are given by (2.37) with fi replaced by On the other hand, 



m = 2^fMJi - (3.8) 

where the factor of 2/x/T, as compared to (2.41), arises from the difference between 
the factors exp(±/i/T). One then obtain the following expressions for the Bose gas 
including both particles and antiparticles^: 

2T 4 

P = —n-, P = 3 P, ( 3 - 9 ) 
7T 

2T 2 

N = (3.10) 
3.1 Relativistic Bose-Einstein condensation 

Let us show that the expression for iVo (2.40) coincides with the thermodynamic 
definition 

N 

N o = y> (3.11) 



5 If we did not neglect indistinguishability of bosons at high temperature, we would obtain, instead 
of (2.46) @, N = ^.Li 3 ( e ^/ T ), where Li v (z) = Y.T^i 2 "/^ is the polylogarithm @, so that, 
for the system including both particles and antiparticles, N = ^r[ij 3 (e' 1 / T ) — Li 3 (e~~ fl / T )]. It then 
follows from the properties of the polylogarithms [Q that, for x = \n\fT « 1, Li^{e x ) — Li^,{e^ x ) = 
\x, so that, we would obtain, instead of (3.10), Nq = /iT 2 /3, which coincides with Haber and 
Weldon's Eq. (1.5). 



13 



where N is the number of bosons in a three-dimensional box of volume V. Since the 
event number density n is, by definition, 



N N 

n ~W)~ VKt 1 

where At is the (average) extent of the ensemble along the g°-axis (as in our discussion 
after (2.44)), one has 

N = nAt. (3.12) 

The equation of motion (2.2) for q° (with M/fiK, the central value of the mass dis- 
tribution, instead of M, which corresponds to a change of scale parameter in the 
expression (2.1) for the generalized Hamiltonian K), 

dq°t Pi 



dr M/ fix ' 

upon averaging over the whole ensemble, reduces to 

At (E) 
T AV M/hk 



(3.13) 



where T&y is the average passage interval in r used in previous consideration. Then, 
in view of (3. 12), (3. 13), one obtains the equation (2.40). 

Since in the particle-antiparticle case, iV re i = N — N, where N and N are the 
numbers of particles and antiparticles, respectively, is a conserved quantity, according 
to the arguments of Haber and Weldon [|K| pointed out in Section I, iVo = N ie \/V is 
also a conserved quantity, so that it makes sense to talk of |iV re i| bosons in a spatial 
box of the volume V. Therefore, in Eq. (3.10) N is a conserved quantity, so that, the 
dependence of fi on temperature is defined by (we assume that N is continuous at 
the phase transition) 

,= ^. (3.14) 

For T above some critical temperature, one can always find \i satisfying (3.3) 
such that the relation (3.14) holds; no such \i can be found for T below the critical 
temperature. The value of the critical temperature is defined by putting = M/2/ik 
in (3.14): 



(3.15) 



For \fj,\ = M/2/ifi-, the width of the mass distribution is zero, in view of (3.4), and 
hence the ensemble approaches a distribution sharply peaked at the mass-shell value 



14 



M/hk- The fluctuations 5m = y (m 2 ) — (m) 2 also vanish. Indeed, as follows from 
(2.24), (2.25) with /i replaced by and (3. 14), (3.15), 




so that, at T = T c , 5m = 0. It follows from (3.16) that for T in the vicinity of T c 
(T > T c ), _ 

M feT / 

5m ~- — J — JT-T c . (3.17) 



We note that Eqs. (2. 45), (2. 46) do not contain explicit dependence on the chemical 
potential, and hence no phase transition is induced. In fact, at lower temperature (or 
small fix) one or the other of the particle or antiparticle distribution dominates, and 



one returns to the case of the high-temperature strongly interacting gas fl39fl . The 
remaining phase transition is the usual low-temperature Bose-Einstein condensation 
discussed in the textbooks. 

One sees, with the help of (3.4), that the expression for n (3.7) can be rewritten 

as 

1 M 

n = — r — Am uT; (3.18) 
2tH fi K 

since at T = T c , Am = 0, it follows that n = at all temperatures below T c . 
Therefore, the behaviour of an ultrarelativistic Bose gas including both particles and 
antiparticles, which is governed by the relation (3.14), can be thought of as a special 
type of Bose-Einstein condensation to a ground state with p^p^ = —(M/hk) 2 (this 
ground state occurs with zero weight in the integral (3.6)). In such a formulation, 
every state with temperature T > T c , given by Eq. (3.6), should be considered as 
an off-shell excitation of the on-shell ground state. At T = T c , all such excitations 
freeze out and the distribution becomes strongly peaked at a definite mass, i.e., the 
system undergoes a phase transition to the on-shell sector. Note that, for n — 0, Eq. 
(3.12) gives At = 0. Then, since (E) ~ T, one obtains from (3.13) that Tav = 
(this relation can be also obtained from (2.44) for Am = 0), which means that all 
the events become particles. 

Since in both (off-shell and on-shell) phases the temperature dependence of pres- 
sure and energy density are the same, and the velocity of sound, c 2 = dp/ dp, is also 
the same, this phase transition is a second order phase transition. As the distribution 
function enters the on-shell phase at T = T c , the underlying off-shell theory describes 
fluctuations around the sharp mean mass. This phenomenon provides a mechanism, 
based on equilibrium statistical mechanics, for understanding how the general off- 
shell theory is constrained to the neighborhood of a sharp universal mass shell for 
each particle type. At temperatures below T c , the results of the theory for the main 
thermodynamic quantities coincide with those of the usual on-shell theories. 
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In order that our considerations be valid, there must hold the relation T c » 
M/hk, which reduces, via (3.15), to 



Nn 



1 / M N 5 



»- — . (3.19) 



7T 2 \\1K 



For M/n K ~ nil, ~ 140 MeV, this inequality yields A" >> 3 • 10 5 MeV 3 . Taking 
A^o ~ 5 • 10 6 MeV 3 , which corresponds to temperature ~ 350 MeV, in view of (2.46), 
one gets T c ~ 550 MeV ~ Am^. 

If hk is very small, it is difficult to satisfy (3.19) and the possibility of such a 
phase transition may disappear. This case corresponds, as noted above, to that of 
strong interactions and is discussed in succeeding paper |39| . 



4 Concluding remarks 

We have considered the ideal relativistic Bose gas within the framework of a manifestly 
covariant relativistic statistical mechanics, taking account of antiparticles. We have 
shown that in such a particle-antiparticle system at some critical temperature T c 
a special type of relativistic Bose-Einstein condensation sets in, which corresponds 
to phase transition from the sector of relativistic mass distributions to a sector in 
which the boson mass distribution peaks at a definite mass. The results which can be 
computed from the latter coincide with those obtained in a high-temperature limit of 
the usual on-shell relativistic theory. 

The relativistic Bose-Einstein condensation in particle-antiparticle system consid- 
ered in the present paper can represent, along with the Galilean limit c — > oo |34{|, a 



possible mechanism of acquiring a given sharp mass by the particles of the system, 
as a phase transition between the corresponding sectors of the theory. Since this 
phase transition can occur at an ultrarelativistic temperature, it might be relevant 
to cosmological models. The relativistic Bose-Einstein condensation considered in 
the present paper may also have properties which could be useful in the study of 



relativistic boson stars ||40|| . These and the other aspects of the theory are now under 
further investigation. 
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